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Abstract. —In the Korean Peninsula, the current distribution of the warm-temperate and 
subtropical vegetation (including many homosporous ferns) is limited to southern coastal areas. 
Paleoecological data suggest that during the Last Glacial Maximum this vegetation retreated to 
glacial refugia putatively located in southern Japan and/or southern China, followed by a post¬ 
glacial recolonization. Two broad scenarios of post-glacial recolonization could be hypothesized: 
extant Korean populations are derived from multiple source populations (i.e., from multiple 
refugia); alternatively, they originate from a single refugium. To test which of these scenarios is 
more likely, we surveyed patterns of genetic diversity in eight ( n = 307) populations of Cyrtomium 
falcatum from southern Korea. We found extremely low levels of allozyme variation within 
populations coupled with high among-population differentiation. These data best support the 
second hypothesis, and indicate that the current genetic diversity may be a consequence of post¬ 
glacial long-distance dispersal events and subsequent founder effects. In addition, restricted gene 
flow among the discontinuous populations of C. falcatum in southern Korea has likely contributed 
to the high degree of among-population genetic differentiation. From a conservation perspective, 
several populations should be targeted for both in situ and ex situ conservation, as C. falcatum 
exhibits a high degree of divergence among populations. 

Key Words. —Dryopteridaceae, CiTtomium , allozymes, conservation, founder effect, glacial refugia, 
homosporous fern, gametophytic selfing. population history, population structure 


Genetic diversity patterns of plant species are shaped by interacting 
historical, biological, ecological, and demographic factors (Nevo et al., 1984; 
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Hamrick and Godt, 1989; Gray, 1996; Duminil et al., 2007). From a historical 
viewpoint, the Quaternary glacial-interglacial oscillations played an increas¬ 
ingly recognized role in shaping the current distribution of plant species and 
thus, their contemporary levels and partitioning of genetic diversity within 
and among populations (Hewitt, 1999, 2000; Hu et al., 2009). For example, 
populations/species that occurred in formerly glaciated regions usually show 
lower levels of genetic diversity than those from unglaciated areas (e.g., glacial 
refugia) through founder effects and bottlenecks as a result of multiple 
stepwise colonization events (Hewitt, 1996; Widmer and Lexer, 2001; Jimenez 
et al., 2010). Thus, the patterns of genetic diversity maintained by the species 
(especially the spatial distribution of genotypes) are often used to infer the 
location of refuges and the post-glacial migration routes from these, and this 
has been particularly fruitful in Europe and North America (Soltis et al., 2006; 
Weiss and Ferrand, 2007, Hu et al., 2009; Hewitt, 2011). 

Ferns have some life-history traits that are strikingly different from seed 
plants and that have potentially significant effects on patterns of population 
genetic variation. First, fern dispersal occurs via haploid spores. Second, their 
gametophytic generation is independent from the maternal sporophytes. 
Third, owing to their small size, fern spores tend to be dispersed much farther 
by wind compared to most seeds (Tryon, 1970, 1972), although this feature is 
analogous to the tiny, dust-like seeds of orchids (Arditti and Ghani, 2000). As 
in many seed plants, however, the majority of propagules fall around the 
immediate vicinity of the parent (Peck et al., 1990). Fourth, since spermato- 
zoids require transport in water, male gamete dispersal distance of ferns tends 
to be very limited (within a few centimeters; Peck et al., 1990). Finally, in 
many homosporous ferns, in the absence of genetic load a single spore could 
produce a sporophyte via intragametophytic selfing (self-fertilization of a 
haploid gametophyte), enabling the successful colonization of new sites 
(Lloyd, 1974; Flinn, 2006; Edgington, 2007; Wubs et al., 2010). Intragameto¬ 
phytic selfing results, in a single generation, in completely homozygous 
sporophytes (Klekowski, 1972; Vogel et al., 1999a), a situation without 
analogue in seed plants. Thus, it has been suggested that genetically 
polymorphic populations could be attributed to the occurrence of multiple 
independent spore dispersal and establishment events over time, whereas 
genetically monomorphic homosporous fern populations are more likely to 
have arisen from single colonists (i.e., single spores; Pryor et al., 2001). 

Habitat specificity and recurrent gene flow of homosporous ferns should be 
regarded as factors determining the degree of population differentiation. Soltis 
et al. (1989) hypothesized that xeric or rock dwelling ferns would exhibit 
higher among-population differentiation than would ferns occurring in mesic 
habitats, due to limited gene flow among isolated rocky habitats. Since then, 

several population-genetics studies have supported this hypothesis (Pryor 
et al., 2001 and references therein). 

In Korea, many ferns are characteristic of the warm-temperate and 
subtropical vegetation, such as Cyrtomium falcatum (L.f.) C. Presl (Drvopter- 
idaceae), a rock dwelling homosporous fern native to southern and eastern 
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Asia, which is taken here as a case study. This vegetation belt currently occurs 
in a narrow zone along the southeastern and southern coast (Yi, 2011). The few 
available pollen and spore records suggest, however, that this warm-temperate 
vegetation was likely absent from the Korean Peninsula during the Last Glacial 
Maximum (LGM, ca. 21,000 yr ago; e.g., Choi, 1998; Chung et al., 2006; Chung, 
2007; Chung et al., 2010; Yi and Kim, 2010), a scenario consistent with 
regional vegetation reconstructions (e.g., Adams and Faure, 1997; Harrison 
et al., 2001; Hope et al., 2004; Prentice et al., 2011). On the southern coast of 
Korea, the onset of the Holocene (ca. 11,000 years BP) and the accompanying 
climatic amelioration were marked by a sudden increase in abundance of ferns 
in the family Polypodiaceae, and an abrupt decline in herbaceous taxa, 
together with the expansion of cool temperate deciduous broad-leaved forests 
(Chung et al., 2010). The first appearance after the LGM of broad-leaved 
evergreen vegetation in the Korean Peninsula was approximately 8,500 years 
ago (Chung, 2011) and somewhat earlier in Jeju Island (ca. 12,000—10,000 yr 
BP; Chung, 2007), which also coincided with a rise of fern spores, indicating 
warmer and more humid conditions (Chung et al., 2010). These paleovegeta¬ 
tion studies suggest that post-glacial colonization either from southern Japan 
(e.g., Kyushu; Fig. 1) or southern China, which harbored glacial refugia for 
warm-temperate vegetation (e.g., Hope et al., 2004; Gotanda and Yasuda, 2008; 
Lopez-Pujol et al., 2011; Qiu et al., 2011), would be much more plausible than 
persistence of warm-temperate and subtropical vegetation in Korean refugia 
during the Pleistocene glaciations. 

Cyrtomium falcatum is an evergreen homosporous fern that usually grows 
on coastal rocky slopes in the warmer parts of south to northeastern Asia 
(India, Vietnam, eastern and southern China, Taiwan, southern Korea, and 
Japan; Iwatsuki, 1992). However, it has become naturalized in many parts of 
the world (including Hawaii, North America, Australia, western and southern 
Europe, Reunion Island, and South Africa) because it escaped from gardens 
(Roux, 2011). The species, 10-60 cm tall, has a short, erect rhizome, and thus, 
it is highly likely that proximally located individuals within populations are 
distinct genets. In southern Korea, C. falcatum usually grows on crevices in 
steep cliffs, rocks, and man-made vertically oriented stone walls near 
seashores, and thus, populations occur discontinuously. Chromosome num¬ 
bers of n = 41 (diploid) or n = 82 (tetraploid) have been reported for C. 
falcatum in Japan (Iwatsuki, 1992). 

Based on the life-history and ecological traits of homosporous ferns, together 
with the information available on the paleoecology of the Korean Peninsula, 
we hypothesize two broad scenarios for the origin of current populations of 
warm-temperate homosporous fern species in southern Korea. If contemporary 
populations were derived from multiple source populations (i.e., from 
multiple glacial refugia), presumably from southern Japan and/or southern 
China, we would expect high levels of within-population genetic variation as 
consequence of the admixture of genetically divergent lineages arriving from 
different refugia (i.e., the ‘melting pot’ effect that has been described for many 
European trees and shrubs; Petit et al., 2003). Regarding among-population 
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Fig. 1 . Upper figure: map of northeastern Asia. Dotted line indicates exposed coastal lines during 
the Last Glacial Maximum (Shi, 2002; Xu et aL , 2010). Lower figure: collection sites of Cyrtomium 
falcatum in southern Korea; CF-1 and CF-2 in Geojae Island. CF-3 and CF-7 (each separated by ca. 
1.5 km) in Oenaro Island: CF~4 in Naenaro Island; CF-5 in Hong Island; CF-6 in Heuksan Island; 
and CF-8 in Haenam-gun (mainland Korea). 


genetic differentiation, either low or high values would be exhibited 
depending on ecological factors (Hamrick and Nason, 1996). Large populations 
that are continuously distributed should exhibit low inter-population 
variation probably due to high recurrent gene flow between adjacent 
populations. In contrast, high genetic divergence would be expected among 
small disjunct populations because of low rates of gene flow between isolated 
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Table 1 . Summary of genetic diversity measures and mean fixation values (F IS ) observed in eight 
populations o: Cyrtomium falcatum. 


Population 

n 

%P 

A 

AR 

(SE) 

H e (SE) 

Fib 

CF-l 

36 

14.3 

1.14 

1.14 

0.044 (0.053) 

0.043 (0.042) 

-0.071 

CF-2 

36 

14.3 

1.14 

1.12 

0.039 (0.042) 

0.046 (0.049) 

0.154 

CF-3 

51 

9.5 

1.10 

1.08 

0.024 (0.036) 

0.019 (0.028) 

— 0.270* 

CF-4 

42 

14.3 

1.14 

1.14 

0.028 (0.026) 

0.027 (0.025) 

0.037 

CF-5 

23 

28.6 

1.29 

1.28 

0.050 (0.041) 

0.071 (0.050) 

0.299* 

CF-6 

18 

4.8 

1.05 

1.05 

0.019 (0.030) 

0.018 (0.030) 

-0.008 

CF-7 

76 

9.5 

1.10 

1.10 

0.044 (0.057) 

0.044 (0.049) 

0.001 

CF-8 

25 

0.0 

1.00 

1.00 

0.000 (0.000) 

0.000 (0.000) 

na 

Average 

38 

11.9 

1.12 

1.11 

0.031 (0.006) 

0.034 (0.008) 

0.030 b 

Pooled samples 

307 

38.1 

1.38 


0.033 (0.018) 

0.069 (0.034) 


Homosporous fems c 


36.1 

1.63 



0.132 



Abbreviations: n, sample size; %P, percentage of polymorphic loci; A, mean number of alleles per 
locus; AR , mean allelic richness based on a minimum sample size of 18 individuals; H ot observed 
heterozygosity; H e , Hardy-Weinberg (H-W) expected heterozygosity or gene diversity; SE, standard 
error; F IS , fixation index within populations; na, not available (because of monomorphism across 
all the loci examined in this population). 

a Denotes significance (P < 0.05) based on permutation (999 replicates) under the null hypothesis 
of F\s ~ 0. 

t 

Non-signiticant Weir and Cockerham (1984) estimate of F\$ over populations. 

Allozyme-based genetic data from Tables 7 and 8 in Li and Haufler (1999). 


populations. Alternatively, if extant populations were established from 
colonizers coming from a single source (i.e., a single refugium), then within- 
population genetic variation would be low because of long-distance dispersal 
associated bottlenecks (e.g., Hewitt, 1996, 2000). Genetic differentiation among 
populations would be high or low depending on rates of contemporary gene 
flow among the Korean populations. To date, these colonization hypotheses 
have not been empirically tested for the warm-temperate and subtropical 
homosporous ferns native to the Korean Peninsula. In this study, we surveyed 
the levels and distribution of allozyme-based genetic diversity in C. falcatum 
to test which of the post-glacial colonization hypotheses is most likely. 
Achieving a better understanding of the genetic structure of this currently rare 
fern in the Korean Peninsula, in addition, will provide guidelines for its 
recovery and management. 


Materials and Methods 

Sample collection.—We collected one leaf segment (pinna) from each 
individual to minimize damage to the plants. A total of 307 individuals were 
sampled from eight populations of C. falcatum from southern Korea, including 
several islands (Fig. 1 and Table 1). All sampled leaf tissue was kept on ice 
until its transportation to the laboratory, where it was stored at 4 C until 
enzyme extraction. 
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Allozyme electrophoresis .—We extracted enzymes by finely cutting leaf 
samples, adding an extraction buffer (Mitton et ah, 1979), and then crushing 
them with a mortar and pestle. Enzyme extracts were absorbed onto 
chromatography wicks and stored in microtiter plates in an ultra-cold 
( —70 G) freezer until analyzed. We conducted electrophoresis on 13% starch 
gels, with three buffer systems. We used a modification (Haufler, 1985) of 
system 6 of Soltis et al. (1983) to resolve alcohol dehydrogenase ( Adh ), 
diaphorase ( Dia-1 , Dia-2), fluorescent esterase ( Fe-1, Fe-2 ), and cathodal 
peroxidase ( Cpx ). We used system 11 of Soltis et al. (1983) to resolve 
glyceraldehyde-3-phosphate dehydrogenase ( G-3-pdh-l, G-3-pdh-2 ), hexoki- 
nase ( Hk-1, Hk-2), isocitrate dehydrogenase ( Idh ), phosphoglucoisomerase 
[Pgi-1, Pgi-2), phosphoglucomutase ( Pgm-1 , Pgm-2, Pgm-3 ), and shikimate 
dehydrogenase ( Skdh ). In addition, we used the morpholine-citrate buffer 
system (pH 6.1) of Clayton and Tretiak (1972) to resolve fructose-1,6- 
diphosphatase ( Fl,6 ) and malate dehydrogenase ( Mdh-1, Mdh-2, Mdh-3). We 
followed stain recipes from Soltis et al. (1983) except for diaphorase (Cheliak 
and Pitel, 1984). We designated putative loci sequentially, with the most 
anodally migrating isozyme designated as 1, the next 2 , and so on. We also 
designated different alleles within each locus sequentially by a, the next b, and 
so on. The observed enzyme banding patterns were consistent with their 

typical subunit structure and subcellular compartmentalization in diploid 
plants (Weeden and Wendel, 1989). 

Data analysis .—We considered a locus to be polymorphic when two or more 
alleles were observed, regardless of their frequencies. We estimated the genetic 
diversity parameters within populations using the programs POPGENE (Yeh 
et al., 1999) and FSTAT (Goudet, 1995): percent polymorphic loci (%P), mean 
number of alleles per locus (A), allelic richness ( AR ) corrected by minimum 
sample size (n = 18 at CF-6, the population with the smallest sample size), 
observed heterozygosity (H 0 ), and Hardy-Weinberg (H-W) expected heterozy¬ 
gosity or Nei’s (1978) gene diversity ( H e ). Except for AR and H n , these 
parameters were also estimated for the total samples as a whole (i.e., at the 
species level). To test for recent decreases in effective population size 
(bottlenecks), we evaluated differences across loci between the H-W H, and 
the equilibrium heterozygosity (H eq ) expected assuming mutation-drift 
equilibrium. H-W H e is not very sensitive to the fate of low frequency alleles, 
whereas H eq is relatively sensitive to population bottlenecks, and declines as a 
result of the loss of such alleles. These differences ( H e — H eq calculated for a 
number of independent loci) were evaluated using a sign test and a Wilcoxon 
sign-rank test under an infinite allele model using the program BOTTLENECK 
(Piry et al., 1999). Since allelic diversity is generally lost more rapidly than H e 
(Nei et al., 1975), recently bottlenecked populations will exhibit an excess of 
H-W H e relative to H^ (Cornuet and Luikart, 1996; Luikart et al., 1998). 

We used the program SPAGeDi (Hardy and Vekemans, 2002) to calculate 
population-level F IS (inbreeding) and its significance level by 999 permuta¬ 
tions under the null hypothesis of F 1S = 0. To measure deviations from H-W 
equilibrium at each polymorphic locus, we calculated averages of Wright’s 
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Table 2. Results of statistical tests for evidence of recent population bottlenecks in Cyrtomium 
falcatum. Numbers reported are P-values of sign and Wilcoxon sign-rank tests conducted using the 
program BOTTLENECK. 


Population 

Sign test 

Wilcoxon sign-rank test 

CF-l 

0.406 

0.188 

CF-2 

0.442 

0.188 

CF-3 

0.659 

0.875 

CF-4 

0.605 

0.938 

CF-5 

0.445 

0.578 

CF-6 

0.519 

0.250 

CF-7 

0.162 

0.125 

CF-8 

na 

na 


(1965) F\s and F$t (deviations from H-W equilibrium of individuals relative to 
their local populations and local populations relative to the total population, 
respectively) following Weir and Cockerham (1984). Using FSTAT, we 
constructed 95% bootstrap confidence intervals (Cl; 999 replicates) around 
means of Fj S and Fst. and considered the observed F\s and Fst to be significant 
when the 95% Cl did not overlap zero. 

To test the overall pattern of genetic structure at the regional scale (i.e., 
isolation-by-distance effects), we conducted a Mantel test (Mantel, 1967) with 
999 replicates, between all pairwise F ST /{1 — Fst) (Fst was calculated 
following Weir and Cockerham, 1984) and the corresponding logarithm 
pairwise geographical distance (Rousset, 1997) under the null hypothesis of 
no spatial genetic structure (regression slope, p = 0). Finally, to determine the 
degree of genetic divergence among populations of C. falcatum, we calculated 
Nei’s (1978) unbiased genetic identity (f) and distance (D) between pairs of 
populations. Using Nei’s D values, we clustered populations into a phenogram 
following unweighted pair-group method using arithmetic averages (UPGMA). 


Results 

Allozyme variation within populations .—Of the 21 putative loci resolved for 
C. falcatum . eight were polymorphic ( Dia-1 , Fl,6, Fe-1, Fe-2, Hk-2, Idh, Pgm-2, 
and Pgm-3). Allozyme variation within populations was extremely low across 
the eight studied populations: mean percentage of polymorphic loci within 
populations (%P) was 11.9, mean number of alleles per locus (A) was 1.12, and 
mean genetic diversity (H e ) was 0.034 (Table 1). Population CF-5 harbored the 
highest allelic richness and genetic diversity [AR = 1.28 and H e = 0.071; 
Table 1), whereas no allozyme variation was found in CF-8 (Table 1). Slightly 
higher levels of genetic diversity were estimated from pooled samples over all 
populations (n = 307): %P = 38.1; A = 1.38; and H e = 0.069 (Table 1). 
Although we did not conduct any bottleneck test on CF-8 because it had no 
allozyme polymorphism, we found no significant indications of recent 
bottlenecks in any of the remaining seven populations (Table 2). 
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Table 3. Allele frequencies for the three loci with the highest degree of population differentiation 
( Fl,6, Pgm-2, and Pgm-3). 


Population 

CF-1 

CF-2 

CF-3 

CF-4 

CF-5 

CF-6 

CF-7 

CF-8 



Allele frequency 


FI. 6 


Pgm-2 



Pgm-3 


a 

b 

a 

h 

a 

0.000 

1.000 

0.194 

0.806 

0.931 

0.000 

1.000 

0.514 

0.486 

0.292 

0.000 

1.000 

0.980 

0.020 

0.765 

0.000 

1.000 

0.083 

0.917 

0.083 

0,565 

0.435 

0.109 

0.891 

0.000 

0.750 

0.250 

1.000 

0.000 

0.000 

0.000 

1.000 

0.704 

0.296 

0.480 

1.000 

0.000 

0.000 

1.000 

0.000 




b 

0.069 

0.708 

0.235 

0.917 

1.000 

1.000 

0.520 

1.000 


Population genetic structure. —Except for CF-3 and CF-5, population-level 
F IS estimates were not significantly different from zero at the 0.05 level 
(Table 1). These results, as well as the non-significant multi-population-level 
Pis (Pis — 0.030; Table 1 and 95% Cl = — 0.254 to 0.401), indicated that 
populations were generally at H-W equilibrium. Deviations from H-W 
expectations due to allele frequency differences between populations were, 
in contrast, significantly high (F ST = 0.543, 95% Cl = 0.218 to 0.703). This 

level of among-population differentiation was largely due to skewed allele 
frequencies at the three loci Fl,6, Pgm-2, and Pgm-3 (Table 3). 

Pairwise Nei’s (1978) / values between populations were high, ranging from 
0.878 (CF-3 vs. CF-8) to 0.997 (CF-2 vs. CF-7) and with a mean of 0.951 ± 0.011 
(SD), which is comparable with the average values reported for other conspecific 
populations of homosporous pteridophytes (average / = 0.911 ± 0.086, N ~ 16; 
Soltis and Soltis, 1989) and of plants overall (average I = 0.950 ± 0.059, N — 
1,572; van der Bank et al., 2001). The apparent discordance between the high 
values of Fst and the high values of /in C. falcatum is simply due to the fact that 
only polymorphic loci are used for the calculation of Fst. whereas both 
monomorphic and polymorphic loci are employed for estimating pairwise Nei’s 
I. The UPGMA phenogram showed that the eight populations were clustered 
largely in accordance with their geographical locations; CF-l/CF-2 and CF-5/CF- 
6/CF-8 (which are located in the eastern and western extremes of southern part of 
Korea, respectively) were clustered separately (Fig. 2). However, we found no 
significant correlation between pairwise genetic differentiation estimates and 
their corresponding between-population logarithm pairwise geographical dis¬ 
tance ((3 = 0.069, R 2 = 0.016, F = 0.277; Fig. 3), indicating that most variation (ca. 
98%) in genetic differentiation was due to factors other than geographic distance. 


Discussion 

Genetic diversity and structure .—Levels of within-population genetic 
diversity are extremely low in C. falcatum (mean population-level estimates; 
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0.065 0.049 0*033 0.016 0.000 

Genetic distance 

Fig. 2. UPGMA phenogram based on Nei’s genetic distances between populations of 
Cyrtomium falcatum. 


%P = 11.9, A = 1.12, H e = 0.034). Although slightly higher values for these 
genetic diversity measures were obtained from pooled samples over all 
populations ( %P = 38.1, A = 1.38, H e = 0.069), A and H e are still lower than 
expected for homosporous ferns (mean species-level estimates; %P = 36.1, A 
= 1.63, H e = 0.132; Li and Haufler, 1999). The low levels of genetic variation in 
the southern Korean populations may be a consequence of post-glacial long¬ 
distance dispersal events and subsequent founder effects (see below for a 
detailed discussion). 

Populations of C. falcatum were generally at H-W equilibrium (multi¬ 
population level Fis — 0.030), a relatively unexpected finding since many 
homosporous ferns have potential for intragametophytic selfing (Klekowski 
and Baker, 1966), which could cause a substantial deviation from H-W 
equilibrium (i.e., a deficit of heterozygotes) within populations. Consistent 
with this expectation, a considerable excess of homozygotes has been found 
within populations of species of Botrychium and Mankvua (Ophioglossaceae), 
which have subterranean gametophytes that obligatelv self-fertilize via 
intragametophytic selfing (McCauley et a/., 1985; Soltis and Soltis, 1986; 
Watano and Sahashi, 1992; Hauk and Haufler, 1999; M. Y. Chung et al. , 2010). 
However, many diploid homosporous ferns exhibit high outbreeding rates (as 
inferred from inbreeding coefficients; Soltis and Soltis, 1989,1992; Ranker and 
Geiger, 2008), and some studies have suggested that they possess mechanisms 
that promote outcrossing in natural populations (Klekowski, 1973; Haufler and 
Gastony, 1978; Haufler and Ranker, 1985; Wubs et al., 2010). Some of these 
mechanisms promote the formation of functionally unisexual gametophytes 
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Pairwise distance (Ln km) 


Fig. 3. Differentiation between populations of Cyrtomium falcatum. Multilocus estimates of 

pairwise differentiation of Fst/(1 ~ F ST ) are plotted against pairwise logarithm (Ln) geographical 

distances in kilometers according to Rousset (1997). There was a non-significant positive 

relationship between pairwise Fs t/(1 F st ) and pairwise Ln geographical distance (r = 0.126, 

P = 0.277). 

through the asynchronous maturation of male and female gametes and the 
control of antheridia initiation by the pheromone antheridiogen produced by 
maturing female gametophytes (Dopp, 1950; Lloyd, 1974; Haufler and Welling, 
1994; Pajaron et al., 1999). This seems to apply for populations of C. falcatum, 
although we do not know which of the above-mentioned mechanisms is 
promoting outcrossing in the Korean populations of this fern. 

Outcrossing plant species usually maintain most of their genetic variation 
within rather than among populations, whereas selfing species show the 
reverse trend (Brown, 1979; Hamrick et al., 1979). Thus, because populations 
of C. falcatum exhibit high inter-population divergence [F s t — 0.543), factors 
other than mating system are probably important in shaping genetic structure 
among populations of C. falcatum. A high degree of genetic differentiation 
among populations has been observed in other homosporous ferns, including 
Adiantum capillus-veneris (Pryor et al., 2001), Asplenium csikii (Vogel et al., 
1999b), Asplenium ruta-muraria (Schneller and Holderegger, 1996), Asple¬ 
nium septentrionale (Holderegger and Schneller, 1994), Asplenium tricho- 
manes subsp. quadrivalens (Suter et al., 2000), Cheilanthes gracillima (Soltis 
et al., 1989), and Sadleria cyatheoides and S. pallida (Ranker et al., 1996). For 
all these cases, patchiness of suitable habitat (which caused restricted gene 
flow) has been proposed as a major driver of population divergence. This 
habitat trait may also account for the high among-population differentiation 
found in C. falcatum in southern Korea, which is primarily due to allele 
frequency differences at three loci (Table 3). For example, six of eight 
populations were monomorphic at Fl,6; of these, the CF-8 population was 
fixed for the allele a, whereas the other five populations were fixed for the 
allele b. At Pgm-2, CF-6 was fixed for the allele a, whereas CF-8 was fixed for 
the alternative allele b. Apart from the low levels of gene flow, genetic drift 
would have been enhanced by small population sizes. Although current 
populations are of moderate size (M. Y. Chung and M. G. Chung, pers. observ.) 
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and we did not find any indications of recent bottlenecks (BOTTLENECK is 
only able to detect those bottlenecks that have occurred within approximate¬ 
ly the past 2N e —4N e generations; Piry et al., 1999), the possibility of older 
bottlenecks should not be dismissed. 

Inference of colonization history of C. falcatum in the southern Korean 
Peninsula. —Since C. falcatum is a member of the warm-temperate and 
subtropical vegetation community, and the Korean populations are at the 
northern edge of the species’ geographic range, one may hypothesize that it 
endured the Quaternary glacial periods at more southerly latitudes. Glacial 
refugia for many elements of warm-temperate and subtropical flora have been 
proposed to occur in southern Kyushu and also in southern Honshu, Japan (see 
Fig. 1; Tsukada, 1984; Hattori, 1985; Matsuoka and Miyoshi, 1998; Aoki et al., 
2004; Gotanda and Yasuda, 2008). Fern spores usually have a high dispersal 
potential; ca. 500 — 800 km and even 3,200 km are suggested as maximum spore- 
dispersal distances (Tryon, 1970,1972). The Tsushima (Korean) Strait was only 
about 10—20 km wide during the LGM and remained relatively narrow until ca. 
14,000-12,000 yr BP (Park et al., 2000; Lee et al., 2008) being therefore easily 
passable. Even during the Holocene the 200 km channel width would have not 
constituted an insurmountable barrier for spore dispersal. Current populations 
of C. falcatum could also have arrived from the more distantly located southern 
regions, as the East China Sea (ECS) was largely exposed until at least 10,000 yr 
BP (Xu et al., 2010). Therefore, migrations from southern China, Taiwan or even 
from some locations offshore in the southern part of the ECS cannot be ruled out 
(see Harrison et al., 2001; Hope et al., 2004; Prentice et al., 2011). 

The low within-population genetic variation for C. falcatum argues against 
the multiple-refugia hypothesis and supports the second hypothesis that 
the contemporary Korean populations of C. falcatum are descendant from 
colonizers from a single glacial refugium, presumably from southern Japan 
and/or southern China. However, we should bear in mind that these two 
proposed scenarios (multiple vs. single source populations) are the two 
extremes of a spectrum of possibility (e.g., some of the extant Korean 
populations could come from a single source, whereas others could originate 
from the admixture of several lineages). Moreover, many factors could have 
altered and/or modeled these “ideal” patterns, such as the number of 
colonization events, the number of propagules arriving at each colonization 
event, and the occurrence of genetic bottlenecks. For example, if population 
sizes have been historically small, random genetic drift since the post-glacial 
colonization events would have lead to low levels of intrapopulation genetic 
diversity even if the populations originated from multiple sources. In this 
latter case, patterns of genetic variation will be hardly distinguishable from 
those expected for species that immigrated from a single refugium. Clearly, 
more species (especially those continuously distributed) should be studied to 
draw firm conclusions about the post-glacial colonization history of warm- 
temperate homosporous fern populations currently occurring in Korea. 

A similar scenario of glacial survival in remote refugia and post-glacial 
recolonization has been proposed for the homosporous fern Dryopteris 
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aemula. Jimenez et al. (2009) reported a total lack of allozyme variation {H T = 

0.000) of this fern in the Iberian Peninsula, which was attributed to founder 

effects during the Holocene expansion. Later, using five microsatellite loci and 

adding one population from the Macaronesian archipelago of Azores, Jimenez 

et al. (2010) found low levels of genetic variation within populations (total 

heterozygosity, fif = 0.447) and a high degree of population genetic 

differentiation (Fgj — 0.520) in D. aemula. Interestingly, the Macaronesian 

population was much more variable than the Iberian ones and, based on these 

findings, the authors suggested that the Azores acted as a glacial refugium from 

which D. aemula spread northeastward and recolonized mainland Europe 

(Jimenez et al., 2010). The role of glacial refugia as sources of plant diversity 

for the post-glacial recolonization in Europe of the Macaronesian Islands has 

been acknowledged in recent years (e.g., Caujape-Castells, 2011; Femandez- 
Palacios et al., 2011; Hutsemekers et al.. 2011). 

In sum, southern Korean populations of C. falcatum exhibit low within-population 
genetic variation, which may be a consequence of post-glacial long-distance dispersal 
events, presumably from a single glacial refugium, and subsequent founder effects. In 
addition, restricted gene flow among the highly specific rock habitats on which C. 
falcatum occurs discontinuously in southern Korea would have contributed to the 
high degree ot among-j >oj ulation genetic differentiation. 

Conservation implications —An understanding of how genetic diversity is 
partitioned within and among populations is critical to design adequate plant 
conservation plans (Godt and Hamrick, 2001; Sun and Wong, 2001). In order to 
preserve a representative sample of the genetic variation, species with high 
population differentiation require the conservation of more populations in 
situ, and also a more extensive population sampling for ex situ conservation. 
Since C. falcatum exhibits a high degree of divergence among populations, a 
relatively large number of populations should be targeted for both in situ and 
ex situ conservation. Using the formula proposed by Ceska et al. (1997), P = 

1 “ (F’st)" (where P is the proportion of genetic variation desired to be 
preserved and n is the number of populations to be sampled/protected), we 
should protect/sample at least four populations in order to conserve > 90% of 
the genetic diversity found in C. falcatum. Considering allelic richness, allele 
frequencies, and the UPGMA phenogram, we suggest that the populations CF- 
1 and CF-7 from one of the clusters and CF-5 and CF-6 from the other cluster 
deserve both in situ preservation and ex situ conservation in southern Korea. 
Thus, these populations should be protected by law (e.g., bv designing plant 
reserves), whereas spores should be collected and deposited in spore storage 
facilities (e.g., by cryoconservation; Ballesteros et al., 2012). 
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